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We consider the nonlinear dynamics of inflaton fragmentation during and after preheating in the 
simplest model of chaotic inflation. While the earlier regime of parametric resonant particle pro- 
duction and the later turbulent regime of interacting fields evolving towards equilibrium are well 
identified and understood, the short intermediate stage of violent nonlinear dynamics remains less 
explored. Lattice simulations of fully nonlinear preheating dynamics show specific features of this 
intermediate stage: occupation numbers of the scalar particles are peaked, scalar fields become sig- 
nificantly non-gaussian and the field dynamics become chaotic and irreversible. Visualization of the 
field dynamics in configuration space reveals that nonlinear interactions generate non-gaussian infla- 
ton inhomogeneities with very fast growing amplitudes. Lhe peaks of the inflaton inhomogeneities 
. coincide with the peaks of the scalar field (s) produced by parametric resonance. When the inflaton 

■ peaks reach their maxima, they stop growing and begin to expand. The subsequent dynamics is 
determined by expansion and superposition of the scalar waves originating from the peaks. Multiple 
wave superposition results in phase mixing and turbulent wave dynamics. Thus, the short interme- 
diate stage is defined by the formation, expansion and collision of bubble-like field inhomogeneities 
associated with the peaks of the original gaussian field. This process is qualitatively similar to 

, the bubble-like inflaton fragmentation that occurs during tachyonic preheating after hybrid or new 

s ! ■ inflation. 

y—i • 

I. INTRODUCTION 

\o . 

■ The origin of matter in the universe from a decaying inflaton held in the process of (p)reheating is a basic feature 
\»0 ' of all realistic inflationary models. If four dimensional effective held theory is sufficient (for specifics of reheating 

in string theory inflation see e.g. 0) this process is described by the non-equilibrium QFT of particle creation and 
thermalization. In chaotic inflation this particle creation typically involves a period of parametric resonance, when 
occupation numbers of Bose particles rapidly become exponentially large . In this case the QFT is well approximated 
by classical field theory, which can be investigated in detail with classical lattice simulations |3|. A full treatment of 
Onj the quantum field theory of the non-linear stages of preheating in controllable models in the limit of high occupation 
ft, ' numbers is in agreement with the classical approximation 

\ The simplest possible inflationary potential contains a massive inflation field V = ^m 2 (p 2 and the simplest preheating 
• model involves a coupling of the inflaton to another field \g 2 4> 2 X 2 - The regime of parametric resonant particle 
production is understood analytically 0. Backreaction of inhomogeneous fluctuations quickly brings the system 
of interacting scalar fields to a strongly non-linear regime characterized by very high occupation numbers. The 
turbulent regime of interacting classical scalar field waves was studied in detail in numerical simulations most 
of which are based on the LATTICEASY code 0, and in the A</> 4 model even analytically with the kinetic theory 
of Kolmogorov-type turbulence |£g. The least understood stage of (p)reheating is the short, violent transition from 
linear preheating to the turbulent stage, which shows anomalies in the momentum space picture, and in the departure 
from gaussian statistics 0. 

Hybrid inflation is another very important class of inflationary models. At first glance preheating in hybrid inflation, 
which contains a symmetry breaking mechanism in the Higgs field sector, has a very different character than in chaotic 
inflation. Preheating in hybrid inflation occurs via tachyonic preheating 0, in which a tachyonic instability of the 
homogeneous modes drives the production of field fluctuations. In hybrid inflation, the decay of the homogeneous fields 
leads to fast non-linear growth of scalar field lumps associated with the peaks of the initial (quantum) fluctuations. 
The lumps then build up, expand and superpose in a random manner to form turbulent, interacting scalar waves 
0,0. A similar picture emeres in new inflation preheating |Toj . Like parametric resonance, tachyonic preheating 
can be interpreted via the reciprocal picture of copious particle production far away from thermal equilibrium, and 
consequent cascades of energy through interacting, excited modes. 

In this paper we investigate in detail the structure of preheating after chaotic inflation in configuration space. We 
find that the intermediate stage between linear preheating and turbulence proceeds via non-linear growth, expansion, 
and superposition of large value field bubbles, similar to what we earlier observed in hybrid inflation. Thus, the 
bubble-like intermediate structure of the non-linear fields is another consistent feature of preheating. However, the 
details of the non- linear dynamics are different: while the bubbles in preheating after hybrid inflation initially occur 
as isolated patches with large spaces in between them, the bubbles that appear in preheating in the model being 
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considered here appear more densely throughout the space and persist for some time as a pattern of standing waves 
before they begin spreading and colliding. 

In section [H] we describe the model we are considering and review the basic nature of preheating in this model, 
focusing on the (well studied) behavior of the fields in momentum space. In section Hill we discuss different diagnostics 
of the system of interacting scalar fields in configuration space and describe the fully non-linear dynamics of inflaton 
fragmentation. In the concluding section we discuss some of the implications of these results, in particular, for the 
generation of gravitational waves from preheating and baryo/leptogenesis from preheating. 



II. PREHEATING AND THERM ALIZATION IN MOMENTUM SPACE 



We consider the potential 



V^mV + ^W (1) 



where (f> is the inflaton and x is another scalar field that is coupled to it. At the end of inflation is a homogeneous 
field which oscillates as <j>(t) w Jt^ mt smm ^ an( i X is a quantum field with eigenfunctions Xfc(i)e _jkx . The temporal 

part Xk(t) obeys an oscillator equation with a periodic frequency u 2 . — (k 2 /a 2 ) + g 2 (j)(t) 2 . The amplitude Xk(t) thus 
undergoes parametric resonance, leading to large occupation numbers of created particles rife. 

Due to the rapid growth of its occupation numbers the field x) can be treated as a classical scalar field. Its 
appearance is described by the realization of the random gaussian field 



X (t,x) = 2 J d 6 k\ Xk (t)\ cos(kx + # fc ) (2) 

i.e. as a superposition of standing waves with random phases 9k and Rayleigh-distributed amplitudes 

P(\Xk(t)\)d\xk(t)\ = e - lxfc i )l ^l^ffll d\xk(t)\, A =< |xfc(0| 2 >■ O ne can use many different quantities to char- 
acterize a random field, such as its variances < a 2 >— J d 3 kk 2n \xk{t)\ 2 7 the spatial density of its peaks of a given 
height, etc. The scale of the peaks and their density depend on the characteristic scale R of the spectrum, which in 
our case is related to the leading resonant momentum fc* ~ yjgrrujjoa 1 / 4 [2]. At the linear stage the phases Ok are 
constant, so that the structure of the random field x stays almost the same. 

Once one field is amplified in this way, other fields that are coupled to it are themselves amplified 5], so within 
a short time of linear preheating (of order dozens of inflaton oscillations) fluctuations of x generate inhomogenoeus 
fluctuations of the field <j>. It is easy to see that fluctuations of <j> will have a non-linear, non-gaussian character. From 
the equation of motion for cj) 

□0 + m 2 4> 2 + g 2 X 2 <i> = , (3) 

we have in Fourier-space 

4> k + 3Hj>k + {(k 2 /a 2 ) + m 2 ) <fe - g 2 Mt) [ d 3 q Xq xl- q , (4) 



where we neglect the term that is third order with respect to fluctuations; <j>o(t) is the background oscillation. The 
solution of this equation with Green's functions |2j shows that <f> fluctuations grow with twice the exponent of x 
fluctuations. It also shows that the fluctuations of </> are non-gaussian. Sometimes this solution is interpreted as 
re-scattering of the particle x<? against the condensate particle 4>o at rest producing Xfc-g an d <frk, x4>o — * X<50- As we 
will sec shortly, this interpretation has significant limitations. 

When the amplitudes of x an d 4> become sufficiently large we have to deal with the fully non-linear problem. The 
field evolution can be well approximated using the classical equation of motion J3J supplemented by another equation 
for x 

□x + .gVx = o. (5) 

Results of simulations of non-linear preheating using the LATTICEEASY program have been reported in many 
earlier papers 0,0- For chaotic inflation, these results were presented in terms of the time evolution of occupation 
numbers nfe(t) or total number density of particles N(t). Figures IT131 show the results of our simulations in these 
familiar terms of Ufe(i) (in combination k 3 ujknk) and N(t), as well as showing the evolution of the field statistics 
(departures from gaussianity) . Here and for the rest of this paper all simulation results are for model Q with 
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m = 10 _6 M p (fixed by CMB normalization) and g 2 — 2.5 x 10~ 7 . The size of the box was L = 10m _1 and the grid 
contained 256 3 points. We also tried other values of g 2 and found qualitatively similar results. 

Figure ^ shows the evolution of the spectra. The spectra show rapid growth of the occupation numbers of both 
fields, with a resonant peak that develops first in the infrared (k ~ fe*) and then moves towards the ultraviolet as a 
result of rescattering. In figure [2] you can see clearly that the occupation number of x initially grows exponentially 
fast due to parametric resonance, followed by even faster growth of the <fi field due to the interaction, in accordance 
with the solution of eq l@J . 
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FIG. 1: Evolution of spectra in the combination k ujkUk of the <f> and \ fields during and immediately after preheating. Bluer 
plots show later spectra. Horizontal axis k is in units of m 
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FIG. 2: Evolution of comoving number density of <f> 
(red, lower plot) and \ (blue, upper plot) in units of 
m 3 . Time is in units of 1/m 
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FIG. 3: Evolution of the ratio </ 2 ) 2 /(/ 4 ), where / 
represents the <f> field (red, solid) or the \ field (blue, 
dashed) and angle brackets represent a spatial aver- 
age, is a measure of gaussianity. This ratio is one for 
a random gaussian field. Time is in units 1/m 



The gaussianity of classical fileds can be measured in different ways. In Figure [3] we show the evolution of the ratio 
(/ 2 ) 2 /(/ 4 ) (kurtosis), which is equal to unity for a gaussian field. During the linear stage of preheating, the field 
fluctuations are a random gaussian field, see eq. (0), reflecting the initial quantum fluctuations that seeded them. 
The inhomogenoeus field </> is generated as a non-gaussian field, in agreement with the solution of eq 10} . When the 
fluctuation amplitude begins to get large, both fields are non-gaussian. During the later turbulent stage both fields 
begin to return to gaussianity. 

Another known feature of preheating is the onset of chaos, when small differences in the initial conditions for the 
fields lead to exponentially divergent solutions: D(t) ~ e At , where D is the distance in phase space between the 
solutions and A is the Lyapunov exponent (see [j| for details). The distance D begins to diverge exponentially exactly 
after the violent transition to the turbulent stage. 

Let us summarize the picture which emerges when we study preheating, turbulence and thermalization in momentum 
space with the occupation numbers n k . There is initial exponential amplification of the field x, peaked around the 
mode fc*. At this stage the x fluctuations form a squeezed state, which is a superposition of standing waves that make 
up a realization of a random gaussian field. Interactions of the two fields lead to very rapid excitation of fluctuations 
of (/), with its energy spectrum also sharply peaked around fc* . To describe generation of cj) inhomogeneities, people 
use the terminology of "re-scattering" of waves. However, there is a short violent stage when occupation numbers 
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have a sharply peaked and rapidly changing spectrum. The field at this stage is non-gaussian, which signals that the 
waves phases are correlated. In some sense, the concept of "particles" is not very useful around that time. In the 
later turbulent stage when nt(t) gradually evolves and gaussianity is restored (due to the loss of phase coherency) 
the picture of rescattering particles becomes proper. As we will see in the next section, gaussianity is not restored for 
some time after the end of preheating. To understand this violent, intermediate stage, however, it is useful to turn to 
the reciprocal picture of field dynamics in configuration space. 



The features in the occupation number spectra rif.(t), namely, sharp time variations, peaks at k ~ fc„, and strong 
non-gaussianity of the fields around the time of transition between preheating and turbulence suggest that we dealing 
with distinct spatial features of the fields in the configuration space. This prompted us to study the dynamics of the 
fields in configuration space. 

The evolution of the fields in configuration space is shown in Figure 0] Each frame shows the spatial profile of the 
fields <f> and \ along a two-dimensional slice of the 3D lattice. A movie that includes many more time frames can be 
found at http://www.science.smith.edu/departments/Physics/fstaff/gfelder/public/bubbles/ Note that 
here (and everywhere in this paper) times are reported in units of 1/m. 

The initial evolution of the fields (t < 100) is characterized by linear growth of fluctuations of x- During this stage 
the fluctuations have the form of a superposition of standing waves with random phases, which make up a random 
gaussian field J3J ■ The eye captures positive and negative peaks that correspond to the peaks of the initial gaussian 
random field x- The peaks in this early stage correspond to the peaks of the initial gaussian random field \. Next 
the oscillations of \ excite oscillations of <\>. The first panel of Figure 0} shows a typical profile near the end of 
this period, just as the oscillations are becoming nonlinear and <p is becoming excited. The amplitude of these <f> 
oscillations grows much faster than the initial x oscillations (see our discussion of eq. (@J) and the oscillations have 
different (and changing) frequencies. The peaks of the (j) oscillations occur in the same places as the peaks of the x 
oscillations, however, as can be seen in the bottom three panels on the left side of Figure We can analyze the 
evolution during this early non- linear period, not using the Fourier- mode description (0J , but instead considering the 
configuration space description equation Q. The interaction term in the equation of motion is approximated during 
this stage by g 2 4>oX 2 > so we have 



where G(s) is the retarded Green's function of the massive scalar field wave equation, s = (t — t') 2 — (x — x') 2 . 
Thus the profile of <p(t, x) is a superposition of the still oscillating homogeneous part plus inhomogeneities induced by 
the Yukawa-type interaction (g 2 <j>o)4>X 2 m the Lagrangian. Since the Yukawa interaction is a short-range interaction 
(defined by the length scale 1/m for a space-like interval G(r) ~ e~ mr ), induced inhomogeneities of (j) appear in the 
vicinity of those in x- 

In the next stage (t ^> 110) the peaks reach their maximum amplitude, comparable to the initial value of the 
homogeneous field </>, and begin to spread. The two lower left panels of Figure shows the peaks expanding and 
colliding. In the panels on the right you can see the standing wave pattern lose coherence as the peaks send out 
ripples that collide and interfere. By t — 124 the fluctuations have spread throughout the lattice, but you can still 
see waves spreading from the original locations of the peaks. Shortly after that time all coherence is lost and the field 
configurations appear to be like random turbulence. 

Figure J5J) shows the distribution of energy density at several points during the evolution. Several points about these 
figures are worth noting. For the parameters we are considering, the gradient energy is subdominant throughout the 
violent rescattering stage and only begins to be significant during the onset of turbulence. Fluctuations in the potential 
and kinetic energy grow in the locations of the bubbles seen in figure |0J, but they are out of phase such that the total 
energy density remains nearly homogeneous. Later, as the bubbles spread and collide, the phase coherence is lost and 
inhomogeneities appear in the total energy density. A movie showing many more time frames of energy density can 
be found at http://www.science.smith.edu/departments/Physics/fstaff/gfelder/public/bubbles/ 

For the first time we calculate the evolution of the statistics of the fields. Figure © shows the field distributions 
(histograms) at various times during the evolution. Initially both fields have gaussian distributions (from their random 
quantum fluctuations), with <fi sharply peaked around <fio and x centered around zero. As non- linear effects become 
important, the statistics of both fields become quite non-trivial. The distribution of the inflaton field becomes at 
times sharply peaked to one side, when the condensate is at an extreme end of its oscillation, and at other times 
bi-modal, with a distinct presence of the homogeneous component plus a significant inhomogeneous component. The 
statistics of the x field also become strongly non-gaussian. Perhaps most surprisingly, the statistics of both fields 



III. INFLATON FRAGMENTATION IN CONFIGURATION SPACE 
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remain non-gaussian for a long time after preheating. At the end of our simulation, at t = 300, the fields were still 
noticeably non-gaussian. During all this time the random phase approximation of interacting scalars is justified only 
to the extent that the distributions approximate gaussians. 



IV. DISCUSSION 



It is convenient to split the process of transition from the homogeneous inflaton condensate to the radiation of 
randomly moving waves into four stages: The first is exponentially rapid growth of small inhomogeneitics that emerge 
from vacuum fluctuations. During this stage fluctuations are linear and the fields are gaussian random fields. The 
second stage is violent backreaction and rescattering of waves, with non-gaussian, nonlinear, nonthermal fluctuations. 
The third stage is Kolmogorov turbulence. During this period the fluctuations are again (nearly) gaussian and energy 
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FIG. 5: Energy density and its components. The horizontal axes represent the same two dimensional slice through the lattice 
as in figure 



gradually cascades towards high momentum modes. Finally, there is thermalization. 

In this paper we used lattice simulations to investigate the second stage of violent field restructuring. We considered 
the model l[T]l. in which preheating occurs through parametric resonance, and examined the evolution of the fields 
in configuration space. The picture that emerged is similar in many ways to what we observed earlier for tachyonic 
preheating in hybrid inflation. For instance in the F-term inflation example where the non-linear potential around 
the bifurcation point is V(a) = —^v<f) 3 + j(j) 4 + j^v 4 , small initial random fluctuations of cf> are amplified by the 
non-linear 3 term. The peaks of the gaussian field thus begin to grow very fast relatively to the surrounding regions 
of 4>. Along those same lines, we found here that oscillations of the x an d <t> fields grow initially at the locations of 
the peaks in the initial x field. For the model Q these growing peaks form a pattern of standing waves that persists 
throughout the linear regime and then begin to spread and overlap as rescattering becomes important. 

The bubble- like structure we see here (and in hybrid inflation) has nothing to do with first-order phase transitions, 
but just with the initial structure of vacuum fluctuations plus non- linear dynamics [l4j . Growth of the individual 
peaks results in the build-up of the scalar field gradients. Subsequent evolution is defined by the expansion of the 
bubbles. The superposition of many almost spherically expanding bubbles leads to decoherence and turbulent motion 
of the scalar waves. 

The main lesson we have learned from this work is that the preheating stage of linear fluctuations and the turbulent 
stage of interacting waves are divided by a short, violent stage of non- linear formation and collision of bubble- like 
large value field regions. 
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FIG. 6: Field distributions. The horizontal axes are field values and the vertical axis represents the frequency of that field 
value on the lattice. The black dots represent the simulation results and the red lines are best-fit Gaussian curves. The fitting 
lines can not be seen in the first frame because they lie directly under the dots. 
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Eventually the fields reach thermal equilibrium characterized only by the temperature. Does that mean that all 
traces of inflaton fragmentation history are erased? There are potential tracers of the non-linear stage of preheat- 
ing related to out-of-equilibrium processes. For instance, people have discussed realizations of baryogenesis at the 
electroweak scale via tachyonic preheating after hybrid inflation [llj, and this process is ultimately related to the 
bubble-like lumps of the Higgs field that form during tachyonic preheating [T^ . Since we now see that fragmentation 
through bubbles can also occur in chaotic inflation, baryogenesis via out-of-equilibrium bubbles can also be extended 
to these models. 

There is another, potentially observable consequence of the non-linear "bubble" stage of inflaton fragmentation. 
Lumps of the scalar fields correspond to large (order of unity) energy density inhomogeneities at the scale of those 
bubbles, R. Collisions of bubbles generate gravitational waves. The fraction of the total energy at the time of 
preheating converted into gravitational waves is significant. We estimate it is of the order of 

^ ~ (RH) 2 , (7) 

Prad 

where 1/H is the Hubble radius. This corresponds to a present-day fraction of energy density Qgw ~ 10 -5 (RH) 2 . 
The way to understand formula (JJJ is the following: The energy converted into gravitational waves from the collision 
of two black holes is of the order of the black hole masses. If the mass of lumps of size R is a fraction / of a black hole 
of the same size, then the fraction of energy converted to gravitational waves from two lumps colliding is /. Scalar 
field lumps at the Hubble scale would form black holes, so in our case / = (RH) 2 . 
The present-day frequency of this gravitational radiation is 

f*^Hz, (8) 

where M = V 1 / 4 is the energy scale of inflation with the potential V. 

For the chaotic inflation model considered in this paper the size of the bubbles is R ~ few/m and at the time they 
begin colliding H ~ m/100, so that the fraction of energy converted into gravitational waves is of the order 10 -3 — 10 -4 . 
This figure is in agreement with the numerical calculations of gravitational wave radiation from preheating after chaotic 
inflation [l3|. 

For chaotic inflation with M at the GUT scale the frequency JSJ is too short and not observable. Gravitational 
waves continue to be generated during the turbulent stage and even during equilibrium due to thermal fluctuations, 
but with a smaller amplitude. It is a subject of further investigation if they can be observed. The most promising 
possibility for observations is, however, generation of gravity waves from low energy hybrid inflation, where / can 
much much smaller. 
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